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Abstract. In this article we present some results concerning natural dissipa- 
tive perturbations of 3d Hamiltonian systems. Given a Hamiltonian system 
x = PdH , and a Casimir function S, we construct a symmetric covariant ten- 
sor g, so that the modified (so-called "metriplectic" ) system x = PdH + gdS 
satisfies the following conditions: dH is a null vector for g, and dS(gdS) < 0. 
Along solutions to a dynamical system of this type, the Hamiltonian function 
H is preserved, while the function S decreases, i.e. S is dissipated by the sys- 
tem. We are motivated by the example of a "relaxing rigid body" by Morrison 
171 in which systems of this type were introduced. 



Introduction 

In his article 11 A Paradigm for Joined Hamiltonian and Dissipative Systems" (see 
7 ), P.J. Morrison introduced a natural geometric formulation of dynamical systems 
that exhibit both conservative and nonconservative characteristics. Historically, 
conservative systems have been modelled geometrically as Hamiltonian systems 
of the form x — PdH where P is a Poisson tensor and H is a smooth function 
(Hamiltonian function). In a Hamiltonian system, the equality dH/dt = dHPdH = 
can be interpreted as conservation of the "energy" H; thus such systems are 
natural models for conservative dynamics (see Nonconservative or dissipative 
systems can also be described geometrically as gradient systems: x = gdS where 
g is a symmetric tensor, S is a smooth function, and dS/dt — dSgdS is typically 
negative definite. Since the function S is not conserved, it may be interpreted as a 
form of energy that is dissipated (or as the negative of "entropy" which is produced 
by the system) (see [5]). 

Morrison's formulation combines these two types of systems into a so-called 
Metriplectic System 

x = PdH + gdS, 

with the additional requirements that H remains a conserved quantity and S con- 
tinues to be dissipated. These requirements can be met if the following conditions 
on H and S are satisfied 

PdS = gdH = 0. 

That is, S is a Casimir function for the Poisson tensor P and dH is a null vector 
for the symmetric tensor g. Morrison applied this formulation to the equations 
for the rigid-body with dissipation and the Vlasov-Poisson equations for plasma 
with collisions. The formulation of dissipative systems as combined Hamiltonian 
and gradient systems has also been studied by 0, 0, ©,[1], and others in 
mathematical physics and control theory. 
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In this article we regard metriplectic systems in K 3 as dissipative perturbations 
of Hamiltonian systems. In the first section we suggest a natural form for the sym- 
metric covariant tensor g that depends only on the differential of the Hamiltonian 
function H, and prove some results about the equilibria of the combined system. In 
the second section we reproduce Morrison's example as a special case, and present 
some other interesting applications. 



1. A class of Metriplectic Systems in R 3 

Let (M, P) be a three-dimensional vector space equipped with a Poisson tensor 
P and standard Euclidean metric h. At each point x we identify T X M and T*M 
with M 3 via the metric h: i>" = hv b . When the context is clear, we will denote the 
dot product on each space, and the pairing between the two spaces (with respect 
to the metric h) by the same symbol, i.e. u" • = u b ■ v b — u" • v b = u l vi. 

For any function H in C°°(M), the vector field £p = PdH defines a Hamiltonian 
system x = £p. Let S € C°°(M) be a Casimir function for this system (PdS = 0). 
We wish to construct a canonical dissipative perturbation 

x = PdH + gdS, 

of this system, so that the following two conditions hold: H = and S < 0, where 
g is a symmetric covariant tensor on M. In other words, we want g and S to satisfy 

(1) gdH = and dS ■ gdS < 0. 

Let (x 1 , x 2 , x 3 ) be local coordinates on M, and write dH = Hidx 1 and dS = Sidx 1 . 
Assume, for now, that each Hi is nonzero. In order for gdH = to hold, the 
following relationships between the coefficients of g must be satisfied. 

g 11 = -g 12 {H 2 /H 1 )-g 13 {Hz/H l ) 

(2) g 22 = -g 2 \H x lH 2 ) - g 23 (H 3 /H 2 ) 
g 33 = -g 3 \H 1 /H 3 )-g 32 (H 2 /H 3 ) 

With these diagonal terms, we can calculate gdS: 

/-(Si/Hi)(ff 12 H 2 + 9 13 H 3 ) + S 2 g 12 + S 3 g 13 \ 
gdS = S l9 21 - (S 2 /H 2 )(g 21 H 1 + g 23 H 3 ) + S 3 g 23 . 
V Siff 31 + S 2 g 32 (S 3 /#3)(5 31 #i + 9 32 H 2 ) J 

In each cotangent space T*M let <Ti(x) denote the i th component of the cross- 
product d x S x d x H and let a be the one-form o = aidx 1 . Then the vector field gdS 
can be locally expressed as 

f-k(9 13 CT2-g 12 <T 3 )\ 

gds= ^(g 12 <y 3 - .g 2 V) . 
U(.9 32 ^ 3 V 2 )/ 
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Therefore, 



1 ',13„ „12„ \ i ( r \2„ „23„ \ , ^3 i 32 „ „31, 



III t±2 tiz 

= „ it (^3-^2 — S2H3) + — {S1H3 — S3H1) + (S2H1 — S1H2) 

Il 2 tl 3 nitLz tl\tl2 

' -(a 2 H l9 32 + a 2 H2 9 13 + <J 2 H 3 g 12 ). 



H1H2H3 

According to the second condition in we must choose coefficients <? y , such that 
this quantity is non-positive, ff we take g^ — H % W for i 7^ j (we have lifted indices 
of dH via the metric h), then we have 

(3) dS ■ gdS = -(erf + of) = - \\a\\ 2 < 0. 

Substituting H % W for g 1 ^ (i 7^ j) into (0) we find that the diagonal terms of g 
should have the form g u — — WHK Thus, we can construct a tensor g that 
satisfies the required conditions: 

(4) //'•' //'//' s't/rir. 

or, invariantly: g — WH (g) WH — l\\WH\\ 2 , where Vif — dH^ and I is the unit 
tensor. The rank of g is zero only at the points for which dH = 0, and since 
gdH = 0, it is never more than two. In fact, these are the only possibilities. 

Lemma 5. // d x H 7^ then the tensor g = VH ® VH — I\\ V if || 2 has rank 2 at 
the point x. 

Proof. Define the following vectors at each point in M: 

Vl = {{),H\-H 2 ), v 2 = {H 3 1 0,-H 1 ), v 3 = (H 2 ,-H\0). 

Observe that at each point in M for which dH 7^ 0, the set J = {vi, V2,v 3 } spans a 
subspace of dimension 2. A simple calculation shows that g(vk)^ — — \\dH || 2 vt for 
each k, so the set J is contained in the image of the homomorphism %g : T*M — > 
TM. Hence, at points for which \\dH \\ 2 7^ 0, the rank of g is 2. □ 

Consider the map (jj<?) b from T*M — > T*M defined on v S T*M by lowering an 
index of %g(v), i.e. \$gf(v) = (gvf . The k th component of {lgf{dS) is 

(gdS) k = (H k W - 5{H % H l )S 3 

= HkH ] Sj — HiWSk 

= (dH ■ dS)H k - [dH ■ dH)S k 

= [dH x (dH x dS)] k . 

Thus, (gdS)^ can be interpreted geometrically as (a multiple of) the component 
of dS which is ^-orthogonal to dH, and so jj<7 at the point x can be interpreted 
as a projection operator from T*M onto the level set of H that passes through 
x (see figure^. The vector field PdH is a Hamiltonian vector field, and so the 
vector Pd x H also lies on this level set. Since gdS is constructed to perturb the 
Hamiltonian system, we would expect gdS and PdH to be independent at most 
points; in fact we can say more that this. 

Lemma 6. At any point in M , the vectors PdH and gdS are h-orthogonal. 
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Figure 1 . Vectors along a solution 



Proof. Since P is skew-symmetric, and since S is a Casimir for P, we have: 

PdH ■ gdS = {PdH) l [dH x (dH x dS)]i 

= (PdH) l [{dH ■ dS)dH - {dH ■ dH)dS] t 

= (dH ■ dS)[PdH ■ dH] - [dH ■ dH)[PdH ■ dS] = 0. 

□ 

Using the degenerate covariant metric g, we define the metriplectic system 
(7) x = PdH + gdS, 

which satisfies H = and S < 0. As for the unperturbed Hamiltonian system, 
trajectories of Q are contained in level sets of H, but they are no longer confined 
to the symplectic leaves defined by P. In fact, by construction, the function S is 
non-increasing along trajectories: dS/dt — dS ■ PdH + dS ■ gdS — — \\o-\\ 2 . Notice 
that this quantity is zero only when a = 0. That is, the derivative of S along 
solutions to (ZJ is zero either when d x S is proportional to d x H (i.e., the level sets 
for H and S are tangent to each other), or when one of dS or dH vanishes. In any 
case, if dS/dt is zero for all t after some time T, then the system has come to "rest" 
in the sense that the trajectory is no longer transverse to the symplectic leaves - 
the dissipation has been turned off. To see this, consider the dissipative part of 10. 
If dH — 0, then g = 0; if dS = 0, then gdS = 0; and if dS is parallel to dH, then 
gdS is again zero. In these "rest states" the system is purely Hamiltonian, and it 
is to these states that (in most cases) the system will tend. 

In order to analyze the behavior of the system as it relaxes to a rest state, as 
well as its behavior once it has relaxed completely (if it does so), we first make 
a simple observation regarding the relation between the equilibria of Q and the 
equilibria of the unperturbed Hamiltonian system. According to lemma ©, we 
know that the vector field PdH + gdS is zero if and only if each component is zero 
individually. In fact, at most points of M, a stronger relation holds. Recall that a 
regular point for P is a point at which the rank of P is maximal. 

Lemma 8. Define the following vector fields: £p = PdH and £ = PdH + gdS. 
Then, at regular points x G M , £p(x) = <-> £(x) = 0. 
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Proof. If d x S = 0, then £p(x) = £(x), so we may assume that d x S ^ 0. Similarly, 
if d x H = 0, then g(x) = and again £p(x) = £(x), so we may further assume that 
d x H^0. 

If £p(x) = 0, then d x H is in the kernel of P which, since x is regular, is spanned 
by the covector d x S. Thus d x S — Xd x H and 

£(x) = Pdi? + gdS* = A.gdff = 0. 

Conversely, if £(x) = 0, then = e2 x <S • £(x) = e^S 1 • gd x S = —d x S x d x H (see © 
above). Since both factors of this product are nonzero, they must be proportional: 
d x H = Xd x S. But then 

£ P (x) = PdH = XPdS = 0. 

□ 

If x is not a regular point for P, then P(x) — 0, and so £p(x) = 0. Such 
points are clearly equilibria for the Hamiltonian system x = PdH, but the vector 
£(x) = gd x S may not be zero, so, in general, the conclusion of the above lemma 
fails to hold in the non-regular case (see EB1 below). Although points of degeneracy 
of P are not typically points of great interest for the Hamiltonian system, they can 
be relevant to the dynamics of the metriplectic system Q. To distinguish regular 
and non-regular points of M we define the following set: TZp = {x|P(x) ^ 0}. 
Equilibrium points that are in TZp will be called regular equilibria. Applying the 
previous lemma to the two systems x = £p and x = £ gives us the following 

Proposition 9. The system x = PdH + gdS and the unperturbed Hamiltonian 
system x = PdH have the same regular equilibria. 

In other words, perturbing a Hamiltonian system in this way does not alter the 
regular equilibria of the system. Since the regular equilibria of a Hamiltonian system 
have a nice geometric description, we can describe fixed points of Q geometrically. 

Proposition 10. The regular equilibria of the system x = PdH + gdS are either 
critical points of H , or are points where the level sets of S and H are tangent to 
each other. 

Proof. The vector field £p = PdH vanishes at a point x £ TZp exactly when either 
d x H = or d x H £ ker(P). In the second case, d x H annihilates every vector 
tangent to the symplectic leaf through x. Thus, d x H is normal (/i-orthogonal) to 
the symplectic leaf through x, which, for x £ Tip, coincides with the level set of S 
through x. □ 

The dissipative system x = gdS has its own set of equilibria, but these are 
not always preserved when the two systems are combined. Points at which the 
dissipative term gdS vanishes also have a geometric description, related to the 
extreme values of the function S on the level sets of H. 

Proposition 11. The vector fields £ = PdH + gdS and £p = PdH coincide either 
at critical points of H , or at critical points of the function S , restricted to a level 
surface H = Hq . 

Proof. The vector field gdS vanishes in one of three cases: dH = 0, dS = 0, or 
dS, dH 7^ with dS £ ker(g). In the last case, since the rank of g is 2, the vectors 
dS and dH must be proportional, i.e. dS — XdH. Thus, £(x) = £p(x) if and only 
if either d x H — 0, or x is a critical point for the constrained function S\p, . □ 
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Comparing these results, we obtain the following relation between the equilibria 
of the component systems x = PdH and x = gdS. 

Proposition 12. If d x S 7^ at a regular point x, then gd x S — if and only if 
Pd x H = 0. 

Proof. If Pd x H = 0, then LemmaEl together with Prop. 03 tells us that gd x S = 0. 
On the other hand, from the proof of Prop. ^Jwe know that if gd x S — 0, with 
d x S ^ 0, then d x H = Xd x S for some A > 0, and so PdH = XPdS = 0. □ 

Corollary 13. If d x S 7^ at a regular point x, then dS/dt = if and only if x is 
an equilibrium of 10. 

Proof. If gd x S = 0, with d x S 7^ 0, then dS/dt — dSgdS = 0. The converse follows 
from Prop. □ 

In particular, if dS 7^ and P 7^ on a level set Hq of H, then the only rest 
states on Ho for the system {7J are equilibrium points. We also have the following 
counterpart to Prop. [5] 

Corollary 14. A regular point x, for which d x S 7^ 0, is an equilibrium point of 
the system x = PdH + gdS if and only if x is an equilibrium of the gradient system 
x = gdS. 

If dS — at a regular point xq in M, then PdH may or may not vanish (see 
examples below). The vector field £ = PdH + gdS reverts to the Hamiltonian one 
PdH at this point, and £ is tangent to the level sets Hq and Sq through xq of both 
H and S. The following result tells us that in this case, dS — along the entire 
trajectory x(t) that passes through x^, and so £ remains Hamiltonian along x{t). 

Proposition 15. Let x(t) be a solution of J7J through the point Xq. If d Xo S = 0, 
then d x ^S = for all time t. 

Proof. Let Sq be the symplectic leaf through xq. The one-form dS is constant along 
Sq since its Lie derivative along any Hamiltonian flow is zero: 

£pdFdS = d(ipdFdS) + ipdF(ddS) = d(dS ■ PdF) = 0, 

for any smooth function F. Since dS = at the point xq, it must remain zero on 
the entire leaf, and so the vector field £ = PdH + gdS reduces to £ = PdH on So. 
Thus, the solution x(t) must be contained in the set So- O 

The symplectic leaf So through a point xq is a submanifold contained in the level 
set So of S, with dimension equal to the rank of P at xq. When xq is a regular 
point, the rank of P is 2, i.e. 1S0 = So- Thus, when d Xo S = 0, with xq regular, the 
vector £ is tangent to So at every point, and any trajectory of J7J that starts on Sq 
must remain there for all time, i.e. the level set So is invariant under the flow of 

The system restricted to Sq is a Hamiltonian system - dissipation is turned 
off. Moreover, since dS = along So, the function S(t) is constant along the flow 
through xo- In this case, the system is in a "rest state"; it is a conservative system 
with constant (typically maximal or minimal) "entropy" . 

2. Examples and Applications 

In this section we discuss several examples of metriplectic systems in M 3 of the 
type 0. We apply the perturbation method developed in the previous section to 
some classical Hamiltonian systems. 
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2.1. Relaxing Rigid Body (Revisited). The example of a "relaxing rigid body" 
by Morrison was the original motivation for our study of metriplectic systems 
of this type. We now reproduce this example via an application of the above 
perturbation method. The Poisson tensor at a point m = (x,y, z) is 

/ z -y\ 
P=\-z x . 

\v o / 

The Hamiltonian and Casimir functions are H = (1/2) (ax 2 + by 2 + cz 2 ) and S = 
(l/2)(x 2 + y 2 + z 2 ). The symmetric tensor g is then 

(—b 2 y 2 — c 2 z 2 abxy acxz \ 

abxy —a 2 x 2 — c 2 z 2 bcyz , 

acxz bcyz —a 2 x 2 — b 2 y 2 1 

which coincides with the operator (up to scalar multiple) defined by Morrison in 
[7j . The vector field PdH + gdS generates a metriplectic system with equations of 
motion given by 

x = (b — c)yz + by(a — b)xy + cz(a — c)xz 
y = (c— a)xz + cz(b — c)yz + ax(b — a)xy 
i = (a — b)xy + ax{c — a)xz + by(c — b)yx 

Notice that these equations can be expressed as rh = PdH + dH x PdH which, 
since PdH = dH x dS, takes the form 

m = dH xdS + dH x (dH x dS). 

The only point at which either dS = or dH — is the origin, which is a 
degenerate point for the system. At any other point, the level sets of H and S 
are ellipsoids and spheres (resp). Along a given level set Ho of H, the equilibrium 
points of the system are points where Hq is tangent to a sphere, i.e. where dH is 
parallel to dS. 

If a 7^ b 7^ c, then this can only occur at the "poles" of Hq where only one of 
x, y, and z is nonzero. At every other point on Hq, since dS/dt = —dS ■ gdS =/= 0, 
S must be strictly decreasing . Thus, the pole with the shortest radius on Hq is 
a stable equilibrium, while the other two poles are unstable (see figure |2J). If two 
of a, b, or c are equal, then two principal radii r*j = rj of Hq are equal, and the 
sphere of the same radius is tangent to Hq along a circle, every point of which is 
an equilibrium. The remaining pole is either stable or unstable depending on its 
length relative to r, (see figure |3J). 




Figure 2. a,b,c distinct 



Figure 3. a = b 
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2.2. Dissipative Oscillators. Here we describe a class of examples in which a one- 
dimensional oscillator x — —x is perturbed by an external function in the following 
way. First, write the system as a two-dimensional Hamiltonian system with the 
standard constant Poisson tensor in R 2 and quadratic Hamiltonian function: 

x = y d_ _ / 1\ (x 

y= -x ° r dt \y) ~ \-l OJ \y 

Introduce an external variable z and construct a new 3d Hamiltonian system with 
Hamiltonian function H = (l/2)(ir 2 + y 2 + z 2 ). 



PdH 



A trajectory through (xq, yo, Za) is either the point (0, 0, zq) on the z-axis (equilib- 
rium), or a horizontal circle at height z$ that lies on the sphere Hq (a level set of 
H) of radius r = \J x 1 . + y 2 + z 2 (see figure 0}. The symmetric tensor g that we 
will use to perturb this system has the form 

(— y 2 — z 2 xy xz 
xy —x 2 — z 2 yz 
xz yz —x 2 — y 2 

Now let S = S(z) be any function of z, and define the metriplectic system m = 
PdH + gdS. In coordinates, 

x = y + xzS' 
(17) y = -x + yzS' 

z = -{x 2 +y 2 )S' 

where S' = dS/dz. The level sets of S are (unions of) horizontal planes, and the 
differential dS is always parallel to the z-axis, except when S'(z) — 0, in which case 
dS = 0. Hence, dS and dH can be parallel only when dH has no horizontal [x or 
y) component, i.e. only at the "north" and "south" poles of a level set of H . 

For a given level set Ho, the poles (0, 0, ±zo), zq > on the z-axis are the only 
equilibria of l|17|l . The stability of each such point depends on the function S. For 
example, if S'(zq) < 0, then any any solution that begins at a point on Hq near 
m = (0, 0, Zq) will flow toward the pole m, and so this point is a stable equilibrium. 

If dS = at some value z = zq, then dS will be zero on the whole plane 
So = {z = zo}. Any trajectory 7 through a point m on So will remain in So for all 
time since the vertical component of its velocity vector will be zero. But 7 must 
also remain on the sphere (level set) Ho containing m, so 7 lies in the intersection 
of So and Hq. This intersection is either a point (when m is a pole of Ho), or a 
circle, in which case 7 is a cycle on Ho at height z = zq. Cycles of this type can be 
either stable or unstable, depending on the function S. 

For example, if we choose S(z) = z 2 , then the equations of motion become 

x = y + xz 2 

y = ~ x + y z2 

2 1 „,2\ 



i = — z(x + y ) 
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Since dS = only at z = 0, the only periodic solutions are cycles at the equators 
of spheres H = k of radius r — y/2k. Clearly z — > along any solution, and 
i 7^ except at z = 0. Thus, these equatorial solutions are attractive cycles for the 
system, and the two poles (0,0, ±r) of each sphere H = k are unstable equilibria 
(see figure QjJ). 




Figure 4. 5 = Figure 5. 5 = z 2 



2.3. Points of Degeneracy. As mentioned above, the system x = Pdff + gdS 
reduces to the dissipative system x — gdS at non-regular points, i.e. points where 
the rank of P is zero. The behavior of the flow at these points varies depending 
on the choice of Hamiltonian function H. Here we present a simple model which 
displays some of the different possibilities in this situation. 

Our model is based on the perturbed oscillator in the previous example, but is 
altered to allow the rank of P to vanish along the (x, z) plane. Specifically, we 
define P as 

/0 y 0\ 

P=\-y o o . 
\ o 0/ 

For a given function H, every point on the plane {y = 0} is a fixed point for 
the Hamiltonian system rh = PdH. For any level set Hq = {H = k} of H, the 
intersection X of Hq and the plane {y — 0} is invariant under the flow of m = PdH . 

We now perturb this system according to the method described above. The 
metric g is given by (|16|) . and the Casimir functions for P are functions that only 
depend on the third coordinate: S = S(z). The metriplectic vector field in this 
case has the following form: 

/ yH 2 + H^S' \ 
£ = PdH + gdS= -yH x + H 2 H 3 S' . 

\ -(H? + H$)S' J 

Which, when y = 0, reduces to the vector field £,\ y= o — gdS, or 

/ H,H 3 \ / H 3 \ ( \ 

C| y =o = S' H 2 H 3 = S'H, + S'H 2 \ H 3 . 

\-(Hf + Hl)J \- Hl ) \-H 2 ) 

The points on the plane {y = 0} are not necessarily fixed by the flow of £, but the 
intersection X will be invariant as long as the vector gdS\ y =o is tangent to X, i.e. 
when the second component of gdS\ y =o is zero. From the expression for £ above, we 
see that X is an invariant set exactly when S' H 2 H 3 \ y=0 = 0. Since S is a function 
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of z only, we have two possibilities: i?2|j/=o = or Hz\ v= q = 0. In the first case, 
we have 

£|y=0 = S Hi 

while in the second case: 

£| y= o = S'Hi ( j + S'H 2 

The following two examples illustrate the behavior of solutions to the equations 
m = PdH + gdS. In the first example, the set X is invariant under the flow of the 
system, while in the second case X only remains invariant under the flow of trivial 
solutions that have initial values on the (x, y) plane. 

Example 1. Let H = (x 2 + y 2 + z 2 )/2 and let S = z 2 /2. Then the vector field £ 
becomes 

i = PdH + gdS= -xy 
\-(x 2 

When y = 0, the tensor P vanishes, and we have 




xz 2 \ 






° 2 1 


= xz 




—xz I 







C|y=o — gdS\ y =o : 

Suppose that m(t) is a solution to rh = £ with m(0) = (xq,0,zo) and let X be the 
intersection of the plane y = with the level set of H that contains m(0). Since 
the y-component of £|y=o is zero, the point m(t) will be in X for all time t. The set 
Hq is a sphere of radius tq — \/ x\ + Zq , so the set X is a great circle on this sphere. 
The points (0, 0, ±ro) are unstable equilibria, as in the examples above, but now 
two new equilibria arise. When the set X meets the points for which dS — 0, the 
vector field £ vanishes. This occurs at the points (±ro,0, 0) on the equator of Hq. 
The point on the positive x-axis is a stable fixed point, while the other is unstable 
(see figure 0). 





Figure 6. X is invariant Figure 7. X is not invariant 

Example 2. Let H = (x 2 + (y - I) 2 + z 2 )/2 and let S = z 2 /2. In this case, the 
vector field £ is 

y(y - !) + xz2 

i = PdH + gdS = ( -xy + (y - \)z 2 
-{x 2 + {y-l) 2 ) Z/ 
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The y-component of gdS\ y =o is zero only if z = 0. If m = (xo,0, zn) is a point 
on the plane {y = 0}, then the flow of £ through m lies on the level set Hq of H 
containing m, but does not remain on the intersection X unless zq — 0. The points 
(a;, 0,0) where X meets the (x,y) plane is an equilibrium point, and these are the 
only points of X which remain invariant under the flow of £ (see figured). 



In this article we examined metriplectic systems of the type {JJ from the point 
of view of perturbations of Hamiltonian systems. We derived a natural form for a 
symmetric tensor g so that the perturbed system x = PdH + gdS dissipates the 
function S while preserving the energy H. We found that a Hamiltonian system 
x = PdH and its metriplectic perturbation have the same regular equilibria, which 
are related to the extreme values of the functions H and S. We also found that, 
for an appropriate choice of 5 1 , the system Q will tend to a Hamiltonian (possibly 
equilibrium) rest state in which the dissipative term vanishes. We presented exam- 
ples of this type of perturbation, including a reproduction of Morrison's example 
of a 'Relaxing Rigid Body'. 

Remark 1. Although our analysis was restricted to three dimensions, it seems 
reasonable that certain aspects of our construction should carry over into higher 
dimensions, including the form of the symmetric tensor g. In fact, metriplectic sys- 
tems have been constructed and discussed in infinite-dimensional settings, examples 
of which can be found in 0] , , and Q . 

Remark 2. The method of perturbation described here can be brought into align- 
ment with the more customary notion of a perturbation by the introduction of a 
continuous parameter that scales the gradient vector field: x = PdH + egdS. It 
would be interesting to study the bifurcations that arise with regard to the stability 
of equilibria in such systems. For a discussion of such a perturbation involving the 
Lorenz system, see [S]. 

Remark 3. This article is a product of the author's doctoral dissertation entitled 
Metriplectic Systems which can be found in the archives of the library at Portland 
State University, or at the link: web.pdx.edu/~djf. 
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